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i—{ ' The shell model is the standard tool for addressing the canonical nuclear many- 

body problem of nonrelativistic nucleons interacting through a static potential. 
^ We discuss several of the uncontrolled approximations that are made in this model 

l , to motivate a different approach, one based on an exact solution of the Bloch- 

Cn 1 Horowitz equation. We argue that the necessary self-consistent solutions of this 

equation can be obtained efficiently by a Green's function expansion based on the 
Lanczos algorithm. The resulting effective theory is carried out for the simplest 
^ ' nuclei, d and "^He, using realistic NN interactions such as the Argonne ulS and 

^Nj , Reid93 potentials, in order to contrast the results with the shell model. We discuss 

QO , the wave function normalization, the evolution of the wave function as the "shell 

f**^ ■ model" space is varied, and the magnetic elastic effective operator. The numerical 

VO ' results show a simple renormalization group behavior that differs from effective 

f^^ ' field theory treatments of the two- and three-body problems. The likely origin of 

^y\ , this scaling is discussed. 

5. 

1 Introduction 

In many text books the shell model (SM) is motivated by the analogy with 
Brueckner's treatment of nuclear matter. While the exact many-body Hamil- 






r^ , tonian operates in an infinite Hilbert space 
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where Tij is the relative nonrelativistic kinetic energy operator and Vij the 
nucleon-nucleon potential, the SM Hamiltonian acts in a restricted space and 
employs a softer "effective" potential, 

HsM = \Y.{T,,+V:/f). (2) 

Motivating Hsm is the notion that the determination of V^^^ might be simpler 
than solving the original A-body problem: the foundation of Brueckner theory 
is that high- momentum contributions to the wave function might be integrated 
out in a rapidly converging series in pnuciear or, equivalently, in the number 



of nuclcons in high-momentum states interacting at one time outside the SM 
space. 

The SM thus represents exphcitly ^ 60% of the wave function that resides 
at long- wavelengths, treating the A-body correlations important to collective 
modes by direct diagonalization. Implicitly the high-momentum components 
are swept into a rather poorly defined "effective interaction," often determined 
empirically. The strength of the SM resides in the first of these two aspects: the 
technology developed for direct diagonalizations is quite remarkable, includ- 
ing recent progress in_Lanczos-based methods EJ, in treatments of light nuclei 
involving many shells B, and in Monte Carlo samplingd'Q. Its weakness is the 
numerous uncontrolled approximations that become apparent when one tries 
to view the shell model as a faithful effective theory (ET). The thesis of this 
talk is that the same numerical strides that have advanced shell model diag- 
onalizations now allow us to remove these uncontrolled approximations. The 
resulting ET has many differences with the shell model and many similarities 
to the effective field theories under discussion at this workshop. 

Among the SM uncontrolled approximations are the following: 

1) Even in lowest order, where only the pairwise interaction of high-momentum 
nucleons is included in H'^^^, the functional form of the resulting effective 
interaction is not as simple as assumed in the SM, 

{\H^"\)sM=(a(3\H""\lS) (3) 

where the Greek symbols label single-particle shell-model states. For exam- 
ple, if the Slater determinants are formed from harmonic oscillator states, the 
two-body matrix elements must carry an additional index labelling the total 
number of quanta in the configuration on which iJ'^-' ' operates a. Thus H'^^-' 
reduces to the shell- model form only when that index is unnecessary, e.g., when 
a lowest-order calculation is restricted to a single shell. Beyond lowest order, 
three-, four-, and higher-body operators are successively added to H^^^ . 

2) Typically W^^^ lacks the symmetries of the original bare 7?, e.g., transla- 
tional invariance and Hermiticity (though the latter is often enforced by hand). 

3) SM wave functions are orthogonal and normed to unity. In ET the effec- 
tive wave functions are naturally defined as the restrictions of the true wave 
functions |^i) to the model space 

\^,)^{SM\^,)\SM) = \^>ff). (4) 

ET 

Thus the norms are less than unity and orthogonality, which holds for the true 
wave functions, is lost when these wave functions are restricted to the model 



space. 

4) Shell model interactions frequently depend on fictitious parameters such as 
"starting energies," introduced to adjust the energy denominator in the two- 
body G-matrix or to account for intermediate-state average energies when the 
two-body G-matrix is iterated to produce some approximation to a higher or- 
der H^ff. 

5) Perhaps most serious, the important issue of effective operators is almost 
never addressed in a meaningful way. In many cases practitioners adopted a 
phenomenological H^^^ which, while successful in producing spectra, provides 
no diagrammatic basis for calculating effective operators or wave function nor- 
malizations. Even in cases where H''-'-^ is derived from some underlying NN 
interaction, the practice is generally to then employ bare operators. In some 
well-studied cases, such as allowed (3 decay in the Ip and 2sld shells, it is 
then recognized that a phenomenological renormalization (e.g., gA —^ 1) of 
operators greatly improves agreement with experiment. But the origin of this 
renormalization and its evolution with momentum transfer q are left unclear. 
The situation is very unsatisfactory and undercuts the shell model as a pre- 
dictive tool. 



2 Self-consistent Bloch-Horowitz Solutions 

We consider the cononical nuclear structure problem of nonrelativistic point 
nuclfions interacting through a realistic NN interaction, such as the Argonne 
ulSO and Reid93El potentials. The question is whether the uncontrolled ap- 
proximations in the shell model can be removed, leaving a more complicated 
but still tractable effective theory. The approach involves three major steps: 

• Formulating a treatment of effective interactions and operators that exploits 
the basic assumption in Brueckner theory — that interactions at high momenta 
can be integrated out in a cluster expansion (essentially an expansion in pa'^, 
where p is the nuclear density and a an interaction range) — but is otherwise 
exact. The convergence of the expansion could then be tested numerically and 
should depend on the operator under study and the momentum transfer. The 
goal would be to distinguish fully converged results from those which require 
higher order calculations. 

• To find numerical tricks for implementing this formulation, demonstrating 
their validity in cases (e.g, A=2,3,4) where the expansion can be carried to all 
orders, so that the answers should then agree with Faddeev and other exact 



methods. 

• To imbed the formulation in a heavier nucleus, where the cluster expansion 
can be carried out only partially. 

There is some reason for optimism that if the first two goals can be 
achieved, the third might yield very accurate results: the Argonne group clus- 
ter variational Monte Carlo effort on -'^^O appeared to yield nearly exact results 
when clusters up to A = 5 were included. In this talk our efforts on the first 
two points will be described. In particular, we will be able to contrast an exact 
ET of the deuteron and "^He with the shell model to illustrate the shortcomings 
of the later: we think the differences are surprising. 

The approach is sketched in Fig. 1. The Hilbert space is divided into a 
long-wavelength "shell model" space, defined by some energy scale Asm, and 
a high-momentum space. One can truncate the latter at some scale Aoo ~ 3 
GeV, characteristic of the cores of realistic potentials, as above this energy, 
excitations make a negligible contribution. That is. 

Hi// {Asm, A^) -^ ^^^^(Asm). (5) 

Aoo large 

All correlations within the "SM" space are included, but the high-momentum 
correlations in the excluded space are limited to n-body, where n is the cluster 
size. Thus the lowest order effective interaction is 

in=2) — "-(n=2fl) ^^ > 

It corresponds to embedding the A-body ladder diagram of Fig. lb between SM 
states: A-2 of the nucleons are spectators, with the remaining pair scattering 
via a two-body ladder. The notation (n=2,0) states that the two-body cluster 
has no explicit dependence on the nuclear density, varying as p^ . The n=3 
A-body ladder of Fig. Ic is similarly 

^(n=3) — -"(ri=2,l) + -"(ri=3,0) V J 

where H'^f^{n = 2, 1) is the two-body part of the three-body ladder 

{al3\HZL2.i)WP') = E («/37l<'£3)l«'/3'7 " «'7/3' + ■••), (8) 

where kp denotes the Fermi level. This decomposition - which is done only 
to emphasize the content of the cluster expansion - illustrates that HJ_„-. 



contains Hf^^^ „s as well as a correction to the two-body interaction that de- 
pends linearly on the density and is obtained by identifying one ingoing SM 
single-particle leg of the three-body ladder with one outgoing leg, summed 
over all occupied states. It also contains a true three-body piece Hf^^^ ^-^ , 
where three SM ingoing single-particle states connect to three distinct outgo- 
ing states after undergoing a series of scatterings outside the SM space. The 
point is simple pedagogy: treatments of successively larger clusters in the high- 
momentum space correct the lowest-order two-body -fff„l2 o~i ^^ adding terms 
proportional to the p, p^ , etc., in the spirit of Brueckner theory. It also adds 
true three-body terms, true four-body terms, etc. Thus an expansion through 



four-body clusters yields H 



order p^, H 
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, the two-body interaction corrected through 



eff 



the three-body interaction through order p, and H,^f^^ 



a density-independent true four-body interaction. 
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Figure 1: Cluster expansion of the eflfectivc interaction. 



The calculation begins with a definition of the "SM" space. The goals of 
handling bound states and of generating an effective interaction that is trans- 



lationally invariant leaves one sensible choice, many-body states constructed 
from harmonic oscillator Slater determinants. To exploit the relative/center-of- 
mass separability of harmonic oscillator Slater determinants, one must separate 
the SM and high-momentum spaces so that all configurations satisfying 

E < KsAfioJ (9) 

are retained in the former. For example, a SM calculation of ^^O with Ksm — 
4-|-Ao, where Aq is the number of quanta in the ^^0 closed shell, would include 
all Ahio configurations, e.g., OpO/i, 2p2/i, and 4p4/i excitations of nucleons from 
the Ip shell into the 2sld shell, Iplh excitations of a Is shell nucleon into 
the 3s2(ilg shell, etc. One can define the projection operator onto the high- 
momentum space by 

Qsm^Q{^smM- (10) 

where h is the oscillator parameter. Thus the included or "SM" space is defined 
by two parameters, Ksm and h. The preservation of translational invariance is 
also important numerically, as it reduces the two-body ladder to an effective 
one-body problem, etc. p. 

The resulting Bloch-Horowitz equationa is then 

H-ff = H + H—^ ttQsmH 

H'"\-^sm) = E\y^sM) \^sm) = (1 - QSA/)|*) (11) 

where |^) is the exact wave function and i/|^) = -B|^). The difficulty posed 
by this equation is the appearance of the unknown energy eigenvalue in the 
equation for H'^^^ . Thus this system must be solved self-consistently. Note 
that there is no explicit reference to the harmonic oscillator in this equation: 
it enters only implicitly through Qsm in distinguishing the long-wavelength 
"SM" space from the remainder of the Hilbert space. 

There is an extensive literature on this and similar equations, often involv- 
a division of H into an unperturbed Hq and a perturbation Hi ^ H — Hq 
13. There are well-known pathologies with this division involving the effects 
of near-by intruder states on the perturbation expansion IIH. Here we explore 
another approach that is nonperturbative and involves, in effect, a computer 
summation of diagrams. The method is based on the Lanczos algorithm and 
offers a remarkably simple solution to the issue of self-consistency. 

In the Lanczos algorithm a basis for representing a Hamiltonian is formed 
recursively in such a way that the resulting Hamiltonian is tridiagonal. Given a 
Hermitian operator H and and an initial normalized vector \vi), the successive 



steps are 



H\vi) = ai\vi) + lii\v2) 

H\V2) = /3i|wi) +a2\v2) + P2\vz) 

H\v'i) = /32|w2) + a'i\v'i) 



- I33\v4) etc. 



(12) 



so that the H takes the form 



H 



/ai /3i 

Pi 0L2 P2 

132 as 

V ; : ; 



(\vi)\ 

V2. 
J \ ■■ J 



(13) 



The remarkable property of this algorithm has to do with truncating the pro- 
cess in Eq. (12) after n steps, where n can be much smaller than the dimension 
of the Hilbert space. The resulting truncated matrix in Eq. (13) then con- 
tains the information needed to reconstruct the exact 2n-l lowest moments 
of H over the eigenspectrum. As extremum eigenvalues are crucial to higher 
moments, one common application of the Lanczos algorithm is in determining 
such eigenvalues and their associated eigenfunctions. Another is to begin with 
the vector \vi) — 0\g.s.) and then use the algorithm to calculate the moments 
of the response of the ground state \g-s.) to the operator O. A small number of 
moments, e.g., ~ 100, often is sufficient to construct a response function with 
a numerical resolution cQinparable to that achieved experimentally. 

A third applicationta is in constructing fully interacting Green's functions. 
One finds 

-^-^|«i) = gi{E)\v,) + g2iE)\v2) + ■■■ (14) 

where the gi{E) are continued fractions that depend on ai,[3i and where E 
appears only as a parameter. For example. 



9i{E) 



1 



E — ai 



JL 



(15) 



E~a2 



=■3-/5; 



It follows that the Bloch-Horowitz equation can be solved self-consistently 
with only a single solution of the effective interactions problem, even in cases 
where multiple bound states are needed. The procedure is: 

• For each relative-coordinate vector in the SM space I7), form the excluded- 
space vector |wi) = Qsa/^It) and the corresponding Lanczos matrix for the 



7 



operator QsmH. Retaining the resulting coefRcients ai,[3i for later use, con- 
struct the Green's function for sonic initial guess for E and then the dot product 
with {i\H to find {i\H''ff{E)\-i). 

• Perform the "SM" calculation to find the desired eigenvalue E' which, in 
general, will be different from the guess E. Using the stored ai,(3i, recalculate 
the Green's function for E' and H^^f{E') then redo the "SM" calculation. Re- 
peat until convergence, i.e., until the input E' in the Green's function equals 
the desired output "SM" eigenvalue. 

• Then proceed to the next desired bound state, e.g., the first excited state, 
and repeat the above step. Note that it is not necessary to repeat the H'^^^ 
calculation. The eigenvalue taken from the "SM" calculation is, of course, that 
of the first excited state. The procedure then generates distinct H'^^f {E')s for 
each desired state. 

The attractiveness of this approach is that the effective interactions part of 
the procedure, which is relatively time consuming as it requires one to perform 
a large-basis Lanczos calculation for each relative-coordinate starting vector 
in the "SM" space, is performed only once. The diagonalization in the model 
space is generally much faster: modern workstations can handle even large- 
dimension shell model calculations (sparse matrices of rf ^ 10^) quickly (~ 30 
minutes). In practice we found that self-consistency is achieved easily: six to 
eight cycles is typical. (More cycles are required for states with small binding 
energies.) Thus it is quite practical to derive the exact H'^f^(E)s for a series 
of bound states. 

Now we discuss the results of applying this procedure to the simplest nuclei, 
d and "^He, carrying the above process to completion (two-body and three-body 
ladders, respectively). The motivation is two-fold: demonstrate the numerical 
procedures we described above, and provide for the first time exact effective 
theory results that can be compared to those of the shell model. 

The harmonic oscillator mode expansion must be sufficient to represent 
both the long-distance tails of bound states and the short-distance "hard core" 
scattering predicted by realistic NN potentials. (The Argonne vl8 ^5*0 po- 
tentials are shown in Fig. 2.) Inclusion of high- momentum states through 
Aoo ^ 50 yields a deuteron binding energy accurate to ^ 60 keV; extending 
this to ^ 140 produces a result accurate to one keV, provided one chooses 
an oscillator parameter that optimizes the convergence. (However, as we will 
discuss in Section 3, there appears to be a simple scaling with Aoo that allows 
one to extrapolate results to Aoo -^ oo. When this is done, the small binding 



energy differences that result at Aqo = 140 for a reasonable choices of b all but 
disappear. We will return to this point in Section 2.) Fig. 3 shows the rate 
of convergence as a function of b and Aqq. The rate of convergence for •^He is 
similar to that for the deuteron: a '^ 60 keV energy accuracy is achieved at a 

Aoo of -- 50. 
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Figure 2: The Argonne vl8 So potentials. 



The binding energies and operator matrix elements for simple systems like 
■^He can, of course, be calculated by other methods, e.g., Faddeev techniques 
or Green's function Monte Carlo. We thus want to stress that the point of the 
following discussion is to do analogous calculations in the context of an effective 
theory, so that we begin to see the shortcomings of conventional techniques 
like the shell model as well as possibilities for overcoming those shortcomings. 
A first test of the techniques outlined above is to solve the Bloch-Horowitz 
equation for some SM-like space to then see if the resulting self-consistent 
energy is, indeed, the correct value. For model spaces of 2, 4, 6, and Shu; in the 
case of the deuteron we obtained a binding energy of -2.224 MeV (using \/2b — 
1.6/ and Aqo = 140). The exact result is -2.2246 MeV. Similar agreement was 
obtained for ■^He. 

More interesting is the evolution of the wave functions, shown in Tables |^ 
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Figure 3: Deuteron ground state energy convergence as a function of A^c for several choices 
of the oscillator parameter b. This b is defined as in the independent-particle SM: the 
corresponding b for relative motion is V2 times the values shown. It is clear that if an 
inappropriate size is chosen for the basis states (e.g., b=1.9), the rate of convergence can be 
greatly slowed. 



and 0. Various ET calculations were done in small model spaces, analogous to 
shell model spaces, consisting of all 2fiLLj configurations, all Ahu! configurations, 
etc. For each such space we then solved the Bloch-Horowitz equation via the 
Lanczos Green's function method described above, iterating the shell model 
calculation until the self-consistent energy is fully converged. The deuteron 
and "^He calculations involve two- and three-body ladder sums in the excluded 
space, yielding sets of two- and three-body "SM" matrix elements of H'^^f for 
the model spaces. The deuteron calculation is rather trivial; for Aco ~ 50 the 
''He BH calculation involves a dense matrix of dimension ^ 10*, still rather 
modest by current SM standards. The matrix is dense because we work in 
relative Jacobi coordinates, .mther than axLJn-scheme, utilizing standard Talmi- 
Brody-Moshinshy methodsEj. (See Ref. EJ for details.) The results in Tabled 
were obtained with approximately 100 Lanczos iterations: it is apparent that 
the convergence is then quite good. The wave functions must be normalized 
according to Eq. (11): this involves calculating unity as an effective operator. 
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Table 1 : ET results for the deuteron ground state wave function calculated with the Argonne 
ril8 potential. The columns on the right correspond to different choices of the ET model 
space, the analog of a SM space. The rows correspond to the resulting amplitudes for the 
designated, selected configurations | n,Lj). The quantities within the parentheses are the 
square of the norm of the effective wave function, e.g., the probability that the deuteron 
resides in the corresponding "SM" space. 



basis state 


amplitude 


Ohuj 


2huj 


Ahuj 


6?iw 


8 huj 


exact 


(65.9%) 


(79.5%) 


(86.1%) 


(91.3%) 


(93.0%) 


(100%) 


1 1, '^i> 


0.81155 


0.81154 


0.81155 


0.81155 


0.81152 


0.81155 


1 2, ^5i) 


0.00000 


-0.31483 


-0.31483 


-0.31483 


-0.31482 


-0.31483 


1 1, 'D^) 


0.00000 


0.19524 


0.19524 


0.19524 


0.19523 


0.19524 


1 3, ^5i) 


0.00000 


0.00000 


0.24945 


0.24945 


0.24944 


0.24945 


1 4, ^^i) 


0.00000 


0.00000 


0.00000 


-0.20851 


-0.20850 


-0.20851 


1 5, ^5i) 


0.00000 


0.00000 


0.00000 


0.00000 


0.12596 


0.12596 



We will return to this point below. 

The tables show the lovely evolution of the wave function in ET, an evolu- 
tion quite unlike that of typical shell model calculations. The wave functions 
obtained in different model spaces agree over overlapping parts of their Hilbert 
spaces. Thus as one proceeds through 2Tiuj, Ahcu, GTicu, ... calculations, the ET 
wave function evolves only by adding new components in the expanded space. 
The normalization of the wave function grows accordingly. Thus, for '^He, the 
OTiui ET calculation contains 0.311 of the full wave function in the effective 
space; the 0-|-2+47ia; result is 0.700. 

This evolution will not arise in the standard SM because the wave function 
normalization is set to unity regardless of the model space. It will also not arise 
for a second reason, illustrated in Table ||. The matrix elements of H'^^^ are 
crucially dependent on the model space: the listed results for •^He show that 
a typical matrix {a\H'^^f\P) changes very rapidly under modest expansions 
of the model space, e.g., from 2fiuj to Ahuj. Yet it is common practice in 
the shell model to expand calculations by simply adding to an existing SM 
Hamiltonian new interactions that will mix in additional shells. We suspect 
the behavior found for ^He is generic in ET calculations: it arises because 
a substantial fraction of the wave function lies near but outside the model 
space (e.g., see Table ^). An expansion of the model space changes the energy 
denominators for coupling to some of these configurations, and moves other 
nearby configurations from the excluded space to the model space. Naively, 
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relative changes in effective interaction matrix elements of unity are expected. 
Now we turn to the question of operators. The standard procedure in 
the SM is to calculate nuclear form factors with bare operators, or perhaps 
with bare operators renormalized according to effective charges determined 
phenomenologically at g^ = 0, using SM wave functions normed to 1. As 
we now have a series of exact effective interactions corresponding to different 
model spaces, we can test the validity of this approach. The results for the 
elastic magnetic form factors for the deuteron and ^He are shown in Figs. 
4 and 5. One sees in each case that by the time one reaches a momentum 
transfer q ~ 2.5//, random numbers are being generated: bare operators used 
in conjunction with exact effective wave functions generate results that differ 
by an order of magnitude, depending on the choice of the model space. This 
is not surprising, of course. If one considers the operation 

e^«VT3|5.s.) (16) 

at momentum transfers > 2kp , where kp is the Fermi momentum, most of the 
resulting amplitude should reside outside the long-wavelength model space, in 
any simple view of the nucleus. That is, the strength resides entirely in the 
effective contributions to the operator. If these components are ignored, the 
results have to be in error. 

Clearly the effective interaction and effective operator have to be treated 
consistently and on the same footing. If O is the bare operator, one finds 

(*/|6|*,) = {-^f^d^^fi^ff) (17) 

where 

d'^ff = {l + HQsM- Itt^MI + t^ ttQsmH) (18) 

and where the effective wave function normalization of \^1 ) and |^f ), 
mentioned earlier, must be determined using the effective operator 1, e.g., 

1 = {■^,\^,) = {^ff\[l + HQsM- ^7^)(1+ ^ ^ QsMH)mff) 

hi — HLJSM hi~ QsMtl 

(19) 
These expressions can be evaluated with the Lanczos Green's function methods 
described earlier. When this is done, all of the effective calculations, regardless 
of the choice of the model space, yield the same result, given by the solid lines 
in Figs. 4 and 5. 

We would argue, based on this example, that many persistent problems 
in nuclear physics — ranging from the renormalization of qa in (3 decay cJ to 
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Table 2: As in Table hi only for ^He. The basis states are now designated somewhat 
schematically as | N,a), where A' is the total number of oscillator quanta and a is an index 
representing all other quantum numbers. 



state 


amplitude 


Ohuj 


2hui 


ihuj 


6huj 


8 hui 


exact 


(31.1%) 


(57.4%) 


(70.0%) 


(79.8%) 


(85.5%) 


(100%) 


|o,i) 


0.55791 


0.55791 


0.55791 


0.55795 


0.55791 


0.55793 


|2,1) 


0.00000 


0.04631 


0.04613 


0.04618 


0.04622 


0.04631 


1 2,2) 


0.00000 


-0.48255 


-0.48237 


-0.48243 


-0.48243 


-0.48257 


|2,3) 


0.00000 


0.00729 


0.00731 


0.00730 


0.00729 


0.00729 


|2,4) 


0.00000 


0.16707 


0.16698 


0.16706 


0.16706 


0.16708 


|2,5) 


0.00000 


0.00566 


0.00564 


0.00565 


0.00565 


0.00566 


|2,6) 


0.00000 


-0.00017 


-0.00017 


-0.00017 


-0.00017 


-0.00017 


|4,1) 


0.00000 


0.00000 


-0.02040 


-0.02042 


-0.02043 


-0.02047 


|4,2) 


0.00000 


0.00000 


0.11267 


0.11274 


0.11275 


0.11289 


|4,3) 


0.00000 


0.00000 


-0.04191 


-0.04199 


-0.04208 


-0.04228 


|4,4) 


0.00000 


0.00000 


0.28967 


0.28978 


0.28978 


0.29001 


|4,5) 


0.00000 


0.00000 


0.01059 


0.01059 


0.01059 


0.01059 


|4,6) 


0.00000 


0.00000 


-0.00213 


-0.00212 


-0.00211 


-0.00210 


|4,7) 


0.00000 


0.00000 


0.00998 


0.01000 


0.01000 


0.01000 


|4,8) 


0.00000 


0.00000 


-0.11319 


-0.11327 


-0.11330 


-0.11335 


|4,9) 


0.00000 


0.00000 


0.08446 


0.08447 


0.08446 


0.08448 


4,10) 


0.00000 


0.00000 


-0.08613 


-0.08626 


-0.08632 


-0.08638 


4,11) 


0.00000 


0.00000 


-0.00210 


-0.00211 


-0.00211 


-0.00211 


4,12) 


0.00000 


0.00000 


-0.00252 


-0.00254 


-0.00256 


-0.00257 


4,13) 


0.00000 


0.00000 


-0.00020 


-0.00020 


-0.00020 


-0.00020 


4,14) 


0.00000 


0.00000 


-0.00010 


-0.00010 


-0.00010 


-0.00010 


4,15) 


0.00000 


0.00000 


-0.00012 


-0.00013 


-0.00012 


-0.00012 



Table 3: Selected BH 3-body effective interaction matrix elements for ^He, in MeV, illus- 
trating the strong dependence on the "SM" space. The Argonne vl8 potential was used. 





2huj 


Ahcj 


6hLU 


8 huj 


(0,1 1 H^ii 1 2,1) 


-4.874 


-3.165 


-0.449 


1.279 


(0,1 H'^ff 2,5) 


-0.897 


-1.590 


-1.893 


-2.208 


(2,1 H^ff 2,2) 


6.548 


-2.534 


-4.144 


-5.060 
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bare 140ha;+ all effective 



bare lOho; 

bare 8hu; . 

bare 6hio 
bare 4hio 
bare Ihco 




q (f ) 



Figure 4: The magnetic elastic form factor for the deuteron calculated with the exact H'^^ ' , 
SM wave functions normalized to unity, and a bare operator are compared to the exact result 
(solid line). When effective operators and the proper wave function normalizations are used, 
all results become identical to the solid line. 



the systeroatic differences between measured and calculated Ml form factors 
at Bates Ej — very likely are due to our naive treatments of operators. It 
should be apparent from the above example that no amount of work on H'^f* 
will help with this problem. What is necessary is a diagrammatic basis for 
generating H'^^^ that can be applied in exactly the same way to evolving 
O'^ff, From this perspective, phenomcnological derivations of H^-^^ by fitting 
binding energies and other static properties of nuclei are not terribly helpful, 
unless one intends to simultaneously find phenomcnological rcnormalizations 
for each desired operator in each q^ range of interest. 

3 Numerical Renormalization Group Speculations 

While the Lanczos Green's function method is reasonably elegant, the calcula- 
tions described above have a "brute force" aspect in that the high-momentum 
ladders must be summed to very large Aqo . In the of the deuteron we carried 
the sums to Aqo = 140 to assure one keV accuracy. The calculation for '^He was 
stopped at Aqo = 52, with a resulting binding energy of -6.842 MeV. As our 
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bare 50ha;+ all effective — 

bare lOho; — 

^>^ bare Shco - 

\\: ~ ^ bare 6hio 

\ bare 4hio 

\ \-- -■-. bare 2ha; — 

Mr ~':^<\ 

'I \\\ .1 \ , m^ 

12 3 4 5 

q (f ') 

Figure 5: As in Fig. 4, only for "^He. 



calculations are otherwise exact, this is an upper bound: a variational principle 
operates as Aqo is increased. The corresponding Green's function Monte Carlo 
result is -6.87 ± 0.03 MeV; Faddeev results for the vl8 potential give -6.895 ± 
0.001 MeV. 

The ^He calculations can (and will) be carried further: the existing ma- 
trices are only of dimension 12000. However, there are alternatives to the 
brute-force approach. One rather simple idea is to evaluate, instead of the 
three-body ladder, the difference of the three-body ladder and the appropriate 
sum over the three corresponding two-body ladders, evolving this difference 
with Aoo- The resulting H'^^^ then would be obtained by adding to this dif- 
ference the two-body ladder. While this is a tautology, our expectation is 
that the difference will converge more quickly than the full 3-body ladder as 
a function of Aqo , since the difference would start with all pair wave functions 
being properly correlated. If one were able to truncate the 3-body difference 
calculation at a lower Aqo than the corresponding 2-body ladders, considerable 
efficiency would be gained. This possibility will soon be studied. 

Another possibility has to do with an observation about the way both 
ground-state energies and matrix elements of H^-^^ evolve with Aqq. Our cal- 
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culations are based on the division of the Hilbert space into two sectors 

< A < Asm Asm < A < Aoo (20) 

While the deuterium calculation was done with a large value of Aqo, we can 
bring it to a much lower value in order to study the evolution with Aqo (see 
Fig. 3). When the numerical results were examined, it was found that they 
scaled with Aoo simply, approximately as an exponential in A^. For example, 
calculations were done for the Argonne wl8 potential and b — 1.7f, with Aoo 
— 46, 50, and 54, to which we fit the three-parameter functional form 

Eg.sXAoo) = -2.175MeV + 3.178MeVe-i-°^5(^-/^°)' (21) 

Thus -2.175 MeV is identified as the binding energy for Aoo — + oo, if this 
numerical renormalization group equation were exact. The deuteron binding 
energy at Aoo ~ 50 is off by more than an MeV. Table shows that 95% 
of this difference can indeed be anticipated by examining the running of the 
results for Aoo ~ 50: the projected valued at Aoo = 140 is correct to ~ 30 
keV. For b = 1.13/, the case converging most quickly in Fig. 3, a similar 
30 keV accuracy is achieved for an extrapolation from still lower Aoo, 34-40. 
Numerically, such an extrapolation can be done quite efficiently: with proper 
coding, it costs relatively little "overhead" to obtain additional results in the 
neighborhood of some starting Aoo- 

The corresponding exercise was done for '^He using results for Aoo = 44, 
48, and 52. A very similar functional form is found 

Eg.s. (Aoo) = -6.906MeV -I- 0.403MeVe~i'"('^~/'**^)' (22) 

While we do not have exact results similar to those in Table H, the resulting 
asymptotic value of -6.906 MeV is in good agreement with the Faddeev value 
of 6.895 ± 0.001 MeV. Such 10 keV accuracy compares well with the Green's 
function Monte Carlo result, with its 30 keV accuracy. 

Now all of this poses an interesting question: our "renomalization group" 
evolution is a nontrivial one, as it involves a truncation in harmonic oscillator 
energies for few-body Slater determinants. Yet it appears this truncation is 
leading to a very simple exponential convergence in results. This can be com- 
pared to effective field theory, where the truncation is made on the range of the 
potential and the convergence is a weaker 1/Aoo- Can one derive and therefore 
understand this attractive aspect of a harmonic oscillator mode expansion for 
bound states? 

To attack this problem, we envision setting Aoo to some intermediate scale 
~ 50 (around 750 MeV) — similar to the numerical examples above — and 
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Table 4: A "numerical renormalization group" study of the running of the deuteron binding 
energy with Aca- The results in the middle column are exact. Those in the right column 
were obtained from Eq. (21), determined from the variation of the Eg,s.{^oc) near Aoo = 50. 
An oscillator parameter is 6 = 1.7/, chosen so that this is a slowly converging case requiring 
significant extrapolation. The potential used is Argonne i)18. 



Aoo 


^calculated (jyjgy) 


^projcctca (j^gy) 


46 


-0.874 


fit 


50 


-1.069 


fit 


54 


-1.247 


fit 


60 


-1.476 


-1.480 


70 


-1.771 


-1.773 


80 


-1.961 


-1.962 


90 


-2.077 


-2.071 


100 


-2.143 


-2.128 


110 


-2.179 


-2.156 


120 


-2.196 


-2.167 


130 


-2.204 


-2.172 


140 


-2.207 


-2.174 



performing, prior to our SM effective theory, a preliminary integration over 
modes above this scale. If we write 



H = Ho + H' and H'^^ = Hq -H H' 



T/Off 



(23) 



where Hq is the harmonic oscillator Hamiltonian and H'^^^ is the appropriate 
effective interaction to use below Aoo, then we find 



H'^^ff ^ H' + H'Gl H'^ff 



(24) 



where G\ is the harmonic oscillator Green's function corresponding to high- 
momentum scattering, e.g.. 



Gljrur,)= J2 



*r(^i)>(*r(^^2)i 



A>Aoc 

G0(fi,f2)- 



E-Ea 

A<Aoo 



(25) 



E-Ea 
In Eq. (24) the high-momentum ladder sum involves terms of the form 

HMiri)H'{n)Gl^ {ri.r2)H'{r2WsM{r2) (26) 
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The harmonic oscillator Green's function has the attractive property that 
it can be decomposed as 



Gl^^GlJr+,r^) 



(27) 



where r_|_ = |ri + f*2| and r_ = |ri — r2|- Now the restriction to large A 
corresponds to short times and thus to short propagation distances. Thus 
propagation lengths should shorten as Aqo is increased. Furthermore the 
overall behavior of the Green's function is governed by the factor e~''+''~ . 
Thus the propagator is most contracted along the smaller of these two co- 
ordinates, r< = min(r^,r-). As the ladder will be evaluated between long- 
wavelength SM states, these observations suggest expanding ^c;^(ri)if'(ri) 
and ff'(r2)^sjv^(r2) in Eq. (26) around r> in a Taylor series in r<. (There is 
a more pedagogicaLdiscussion of this technique and its relevance to effective 
field theory in Ref. Q.) 



^ x=high 

> 

^ momentum 




Figure 6: Contraction of a ladder diagram involving only very high momentum excluded 
states to a local operator. 

After a bit of algebra, the net result is the contraction of the ladder sum 
to a local "rainbow" diagram, as illustrated in Fig. 6. The contracted Green's 
function has a form familiar to most effective field theorists, 



Ga„ 



9°ir>) + ^r,g\r>)Wr 



(28) 



where this is to be inserted between wave functions evaluated at r>. The 
function g°(?'>) is obtained by integrating Ga^^ (f + , ?"- ) over r<, while g^{ry) 
involves a second moment in the small coordinate r< , etc. Effectively this sums 
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the ladder according to a local density approximation (g*^), in an approximation 
that also takes into account the gradient of the local density (5^), etc. 

One of the most attractive features of the harmonkipscillator is that the full 
Green's function was recently derived in closed formlla, so that g°(?'>) can be 
calculated by integrating the last expression in Eq. (25): an infinite sum over 
high-momentum modes can be replaced by a finite sum over low-momentum 
modes. The full Green's function is shown in Fig. 7, and the integration that 
produces 5° is depicted in Fig. 8. The resulting g° for Aoo = 30 and 50 are 
shown in Figs. 9 and 10 (dashed lines) along with a similar contraction of 
the full Green's function G". One finds that the dashed lines - the residual 
contribution of very high momentum excitations - rapidly oscillates, with a 
frequency very close to A^. Thus a harmonic oscillator expansion converges 
not because every point in coordinate space is well represented by the included 
states, but because integration of a product of such oscillations with a smooth 
Hamiltonian yields a small remainder. Note that increasing Aqo also extends 
the range in r> spanned by the harmonic oscillator basis functions: extended 
wave functions require more work. 

We are now in the process of implementing this local approximation to 
see how far the integration scale Aoo can be lowered without loss of accuracy. 
Unlike effective field theory treatments, our Taylor series expansion involves a 
known potential and can be carried out term by term until exhaustion sets in. 
The results may provide some estimate of the scales that effective field theories 
need to reach to become very accurate. 

We will close with two speculations. If one folds the dashed-line functions 
of Figs. 9 and 10 with the potential of Fig. 2, the result is not dissimilar to 
a jo(Ar>). Such a function, integrated with a product of harmonic oscillator 
wave functions, yields a factor e"°^ . Thus we are hopeful that this line of 
work may indeed explain the numerical scaling we have previously discussed. 
Our goal is to add to Table H a fourth column in which our phenomenological 
scaling function is replaced by a derived one that yields improved results: the 
"break" seen in Table Hat Aoo ~ 100 (a scale characteristic of the Argonne vl8 
hard core) suggests there is a bit of physics missing from our phenomenological 
function. Second, this suggests that the rapid convergence of the harmonic 
oscillator mode expansion, as compared to effective field theory scaling of 1/Aoo 
which results from absorbing short-range contributions to the NN potential 
into contact terms, may be because it offers a better compromise between 
coordinate and momentum space. Hopefully we will be able to be more specific 
soon. 
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Figure 7: The full (summed over all modes) harmonic oscillator Green's function. For parity 
conserving interactions, only the even or odd projections are needed. The even projection is 

shown. 

/ 



r_ = r> / 






integration: 
G(r+,r_) ^go(r>) 



Figure 8: The integration path in the {r^,r-) plane that is used to contract G(r^,r—) to 
So(''>)- 
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Figure 9: The dashed line is the (?o(^>) that results when excitations A > 30 are summed. 
The solid line is the corresponding result for the full Green's function. It is the rapid 
oscillation of the dashed residual — rather than the accuracy of the harmonic oscillator 
expansion at all values of r> - that is responsible for the convergence. The oscillations are 
very regular and quite close in frequency to A ^ 30. 



4 Outlook and Future Plans 

The ability to efficiently solve "^He as an effective theory in a SM-like space is, 
in itself, quite significant: this means it is relatively straightforward to execute 
a faithful BH treatment of heavier nuclei through order p in both the effective 
interaction and effective operators. The numerical effort is comparable. 

The example of "^He also suggested that many of the uncontrolled approxi- 
mations made in the SM cannot be justified. Recently sophisticated numerical 
advances have been made in SM calculations. We have argued that it may 
be equally important to train this numerical power on the shakey foundations 
of that model. A specific example is our use of the Lanczos algorithm - a 
favorite among SM practitioners - to generate the Green's function, leading to 
a tractable strategy for solving the BH equation self-consistcntly. 

We believe that, if the program in Section 3 has modest success, we should 
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Figure 10: As in Fig. 9, only for A > 50. 



be able to handle ''He with an accuracy similar to that achieved already for 
■^He. That would imply that it is now possible to handle effective interactions 
and operators in finite nuclei through order p^. At that point the exciting 
challenge will be to determine where convergence in the cluster expansion is 
being approached. 

We would like to think that the work reported here will form a bridge 
between the ideas growing out of effective field theories, and the successful 
phenomenology we have achieved in traditional nuclear physics (particularly 
in modeling the NN potential and in solving few-body systems with nearly 
exact methods). It may also stimulate effective field theory, which has not yet 
had an impact on nuclear problems beyond NN and three-body systems. It is 
not inconceivable that a marriage of EFT - which strives to handle few-body 
problems rigorously, including relativity, fundamental few-body forces, etc. — 
with some kind of cluster expansion of the Brueckner type might someday yield 
a controlled theory of finite nuclei. Furthermore, the attractive exponential 
convergence of our harmonic oscillator mode expansion is also EFT "food for 
thought." 
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